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Since the renewal of this grant (NAGW-9) in January 1991, the research has pro- 
gressed smoothly according to the plan in the original proposal with some modifications 
(see Appendix A). The significant achievements during this period are described in the 
following: 

I. Work completed (published and accepted for publication) 

1. S. T. Wu and Ai-Hwa Wang, ’’Numerical Simulation of Extended Corona”, 1991, Adv. 
Space Res., Vol. 11, No. 1, 187. 

2. Tyan Teh and S. T. Wu, ” Model Calculations of the Rising Motion of a Prominence 
Loop”, 1991, Solar Physics, 132, 335. 

3. S. T. Wu, M. T. Song, P. C. H. Marten, and M. Dryer, ’’Shear Induced Instability and 
Arch Filament Eruption: A Magnetohydrodynamic (MHD) Numerical Simulation”, 
1991, Solar Physics, (in press). 

4. S. T. Wu, M. T. Song and E. Tandberg-Hanssen, "A Numerical Simulation of At- 
mospheric Responses Due to Emerging Flux From Sub-Photospheric Layers”, Bull. 
American Astronomical Society, (Abstract) Solar Physics Division Meeting, Huntsville, 
AL, April 1991. 

5. M. T. Sun, S. T. Wu, F. S. Weng, and M. T. Song, "Further Numerical Tests for the 
Progressive Extension Method (PEM) for Extrapolation of Solar Force- Free Mag- 
netic Field”, Bull. American Astronomical Society, (Abstract) Solar Physics Division 
Meeting, Huntsville, AL, April 1991. 

A copy of these reprints and abstracts are attached. 

II. Work In Progress 

The work currently in progress is summarized in the following: 

II.l. Critical shear study 

It is well-known from both observation and theory that photospheric shear motion 
has played a key role in physical processes of the energy build-up and release for the solar 
flare. In order to further our understanding of the occurrence of solar flares it is necessary 
to investigate the triggering mechanism. One popular scenario for the onset of eruptive 
solar flares is that in response to photospheric shear motions the magnetic field evolves 
slowly through a series of magnetohydrodynamic-equilibria until a threshold is reached 
where magnetohydrodvnamic (MHD) non-equilibrium sets in. Thus, a magnetic eruption 
occurs, causing the solar flare. 

To substantiate this claim we have employed our newly developed three- dimensional, 
time-dependent MHD code with gravity to simulate the evolution of the coronal field. 
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We use plasma j3 = 0.1 to closely approximate the condition in the actual corona. Some 
preliminary results are shown in Figures 1 and 2. Figure 1 shows the time sequence of the 
eyolving^rnagnetic field configuration in three-dimensions due to the photospheric motion. 
Figure 2 shows the density enhancement and magnetic field line projected on the x — z 
plane. In this presentation (Figure 2), we observe the following interesting features: 

(i) MHD waves (fast and slow modes) were generated due to the photospheric shear 
motions. Fast-mode waves propagate upward and slow mode waves propagate along 
the field fines. 

(ii) A current was also generated and moves upward slowly. 

Presently, we have begun to analyze these results in detail to reveal the physical 
processes in these data. 

II. 2. Improved PEM for the extrapolation of the nonlinear force-free field 

Since the publication of the first version of our extrapolation scheme (Wu et al. 1990), 
called '’Progressive Extension Method (PEM) 7 ’, further improvement was made on the 
regularization techniques. This new regularization technique can be summarized as follows: 
the expression for averaging given by Eqs. (3.5) and (3.6) of the paper by Wu et al. (see 
Appendix) is modified as follows: 

The regularized solution B X J becomes 


^ = (1-7)^+T^j (1) 

where B UJ is the same as Eq. (3.5) and 


7 (~) = l'o€xp[a{z/ z max )} ( 2 ) 

with 7 0 being arbitrary and which will be determined by trial and error method until 
the solution converges and 


ci — ln( 7 ( ,, m a x ) / 7 a ) ( *-> ) 


Bv setting our goals on 7 and 70 , a is determined. 

Using this modified technqiue. we perform a number of tests by using a rather com- 
plicated analytical nonlinear force-free solution I Low and Lou. 1990) as the input to the 
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Figure 2 The evolution of density enhancement and magnetic field lines projected on the x - z 
plane. 
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numerical model. Figure 3 shows the field lines 1 view from the top generated by (a) analyt- 
ical solutions (Low and Lou, 1990) and (b) numerical solution obtained from extrapolation. 
Figure 4 shows three-dimensional field lines corresponding to the solutions given in Figure 
3. By observing these results, it clearly demonstrates that the PEM is a reasonable method 
to obtain nonlinear force-free fields by using the vector magnetograph data. The accuracy 
of the extrapolation can be shown as a dependence of height. It is recognized from these 
results that the height of the extrapolation is about one tenth of the horizontal boundary 
which is typically about 30,000 - 50,000 km. 

It is understood that the analytical solution is noise free. In order to show that the 
PEM is a viable method to analyze observed vectoral magnetograph data we have chosen 
the 1989 March 10 flare data for our tests. Initial results show some very encouraging 
signs. Since there is no analytical solution to compare with the PEM’s results we shall use 
the potential field model for benchmark testing. Namely, we use the 1989 March 10 data 
to perform potential field calculations using both PEM and Schmidt’s (1964) numerical 
code. These results are shown in Figures 5 and 6. In examining these results, we may 
conclude that the results obtained by PEM and Schmidt are identical. This definitely 
shows the accuracy of the PEM. Since the credibility of the PEM is being established, we 
shall compute the nonlinear force-free field configuration for the 1989 March 19 data. These 
are shown in Figure 7. In comparing these nonlinear force-free results with potential field 
results obtained by the same data, we immediately recognize different physics has been 
revealed by these two different models. For example, there are significant shears exhibited 
in the flare production region in Figure 7, but, as expected, shear cannot be accommodated 
in Figures 5 and 6 as the nature of the potential field model. 

In summary, we feel that it is time to apply the PEM for data analysis with caution. 
Certainly, improvements and modifications for better accuracy will be made during these 
studies. 

II. 3. Study of emerging flux from sub-photospheric layer 

It is believed that the additional energy stored in the solar atmosphere could have 
resulted due to the emerging flux from the sub-photospheric layer. It is difficult to model 
these physics because the characteristics of this part, of the solar atmosphere have several 
orders of magnitude variations in density, pressure and plasma 3 within a very thin layer 
( ~- 2000Ln? ). hence little progress was made. In the past few months, we have been able to 
construct such a numerical model using the Lagrangain method. The Lagrangain method 
is chosen because the usual used Eulerian method cannot accommodate large gradient 
numerical computations. 

Some results are obtained, the steady state density, temperature, pressure and plasma 
beta \ 3 ) distributions as a function of heights ii.e from sun's surface up to 2000 km) 
are shown in Figure S. After we obtained this steady '•'ate structured atmosphere, we 
introduced an emerging flux at the lower boundary and unputed the evolutionary state 
of plasma properties. These results are shown m Figures 9. 10 and IT. These figures 
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3 Comparison between the (a) analytical solution (Low and Lou, 1990) and (b) numeri- 
cal solution obtained by PEM for projected magnetic field liens on x — y plane viewed 
from the top. 



(b) 


Figure 


Comparison between (a) analytical solution (Low and Lou. 1990) and (b) numerical 
solution obtained by PEM for three-dimensional representation of magnetic field lines. 
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represent the temperature, density, pressure, magnetic field lines, velocity and Lagrangain 
grid distribution at times, 50s, 500s, and 1050 s after introduction of the magnetic flux 
respectively. Presently, we are analyzing these results in detail. At first glance at these 
results, we could summarize our present understanding as follows: 

a. The emerging flux leads to the formation of a current sheet at the interface of the old 
and new magnetic field. 

b. This current sheet propagates upward. 

c. The induced plasma flow is oscillating at the Brunt-Vaisala frequency and has a period 
of 240 sec. 

d. Also, we show that there is no oscillation when the gravity is ignored. 

e. The maximum downward flow occurs in the neighborhood of legs of the magnetic loop. 



Figure 8 The distribution ot density, temperature, pressure and plasma^ as a function of height 
in the region of photosphere and chromosphere. 
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III. Tasks for the Third Year Grant (January 1, 1992 - December 31, 1992) 


In view of the present progress, we shall continue the present course and fine tuning 
the present studies according to our schedule presented in the original proposal with the 
modification presented in January 24, 1991. Specific tasks for the coming year (Jan. 1 - 
Dec. 31, 1992) are: 

III. 1 . Radiation and thermal conduction effects on nonlinear evolution of critical equilibria 

This is the task we deleted in FY 1991 because of insufficient funding. Since the suc- 
cessful study described in Section II. 1 in which the radiation and thermal conduction were 
ignored. We felt that if we wish to understand the MHD ” non-equilibria” quantitatively, 
it is imperative to study the higher order effects on the energy transport such as radiation 
and thermal conduction. Thus, we propose to include this in our tasks of FY 1992. 

111.2. Quantitative analysis of the 1989 March flares 

As we discussed in Section II. 2 PEM is a reasonable tool to extrapolate the realistic 
solar magnetic field from the photospheric level upward to the lower corona. With further 
development of PEM, we expect that the PEM will be able to obtain the field configuration 
up to 10 solar radii. Presently, we have obtained magnetic field data from Big Bear 
Observatory and Beijing Observatory during the flare period March 8-15, 1989. We shall 
use the PEM to extrapolate solar magnetic fields (see example shown in Figure 7). We 
expect to obtain the quantitative description of the three-dimensional flare characteristics, 
such as, the pre-flare and post-flare state field configurations, critical currents intensity in 
relation to flare onset, and amount of energy stored and released. 

111. 3. Study the energy and momentum transport from the sub- photospheric layer 

Using the newly developed model described in II. 3 we will study the energy and 
momentum transport from the sub-photosphere upward to the corona to understand the 
physics of corona heating and solar wind acceleration. 

111. 4. Assessment of the three year study 

We shall summarize our new findings, success, failures and the present understanding 
of the solar eruptive phenomena. 

References 
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An Affirmative Action/Eauai ODDorturmy insmution 


Jan. 24, 1991 


Tasks to be performed during the period 1 Jan 1991 - 31 Dec 1991 for 
the Grant NAGW-9. 

1. Critical Shear Study 

(a) Employ our three-dimensional magnetohydrodynamic model 
to study the evolution of critical equilibria for a 
non-linear force-free field. 

(b) Determine the critical values of the physical 
parameters for the solar flare or other eruptive 
phenomena based on the results obtained from (a) 
together with the analytical model described in 
Section II. 1 of the proposal submitted June 8, 1990. 

2. Numerical Tests of the Progressive Extension Method (PEM) 
for the Extrapolation of Nonlinear Force-Free Field. 

We plan to perform further numerical tests to establish the 
limits and accuracy of PEM developed by the principal 
investigator (Wu et al. 1990, Ap . J., 362 . 698-708) for the 
extrapolation of nonlinear force-free magnetic fields from 
the solar vector magnetograph. 

Additional funding of $7,000 is requested to provide support of 
a student for 4 1/2 months to prepare 1989 March Flare data for 
analyses. We expect to employ our nonlinear force- free model to 
deduce the flare characteristics from 1989 March flare data which 
will give us the pre- flare field configuration, post-flare field 
configuration, critical current intensity in relation to flare 
onset, amount of energy released and amount of energy stored. 
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NUMERICAL SIMULATION OF EXTENDED 
CORONA 

S. T. Wu and Ai-Hwa Wang 

Center for Space Plasma and Aeronomic Research and Department 
of Mechanical Engineering , The University of Alabama in Huntsville , 
Huntsville , AL 35899 , U.S.A. 


ABSTRACT 

A three-dimensional, time-dependent magneto hydrodynamic (MHD) model is presented for the study of 
coronal dynamics. The model, written in spherical coordinates, extends from the solar surface (1 f?,, where 
1 R s = 6.95 x 10 5 km) to 15 R 3 . This model was developed with two major issues in mind, namely 
for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to 
achieve these objectives we have employed two different numerical techniques to seek solutions for these two 
different, but related, problems; steady state structures and evolutionary structures. These two numerical 
techniques are: (i) relaxation technique for steady state structures: and (ii) FICE ( Full- Implicit- Continuous- 
Eulerian) technique for evolutionary structures. 

To illustrate this model, we present numerical results for examples of both the steady state and evolutionary 
structure of the corona. These results show the additional physical features which cannot be shown by a 
two-dimensional model. Finally, on the basis of the exploratory calculation, we outline some interesting 
physical features which can be considered for the observing programs of future space missions such as 
SOHO, OSL, CORONAS, etc. 

I. INTRODUCTION 

Since the Skylab-ATM experiments in the seventies, we have recognized that the corona is in a transient 
state in contrast to the previous understanding whereby the corona is always in a quiet orderly state 
(Billings, 1966). It is also further realized that the relationship between the flare and the coronal mass 
ejection is not as consistently intimate a s originally thought ( Hildner et al. 1976). In order to understand 
the physics of this fascinating phenomena of so-called ” coronal transients”, a number of theoretical models 
has been presented in the literature (Hundhausen et al. 1984). All of these theoretical models are based on 
magnetohydrodynamic theory. The methodology used to treat these theoretical models could be classified 
into two categories: (i) analytical methods and (ii) numerical methods. Those models treated by analytical 
methods have to conform to certain strict conditions in which a full description of nonlinear dynamical 
behavior is difficult to achieve; nevertheless, the solutions are exact. On the other hand, the models treated 
by numerical methods could obtain globed descriptions of nonlinear dynamics, but these descriptions are 
not unambiguous and may mislead the physical interpretations. A further limitation to these two categories 
is the fact that all of these models are confined to a two-dimensional geometry. Thus, it is inevitable that 
some arguments in the interpretation of observations have taken place. 

In this paper, we present a newly developed three -dimensional* time-dependent, magnetohydrodynamic 
model for an extended corona. We will suggest that this model could be used to understand the physical 
processes from the comparison of this model’s results with observational data. The theoretical description 
of the model presents the basis for the addition of dissipative mathematical terms that could be used 
to understand additional physical processes from specific observational data. The theoretical description 
of the model are included in Section II. The numerical results are presented in Section III. Finally, the 
concluding remarks are included in Section IV. 

II. ANALYSES 

Mathematical Model 

In this study, we have assumed that the solar atmosphere behaves as a single fluid with negligible dissipative 

effects. With these assumptions, the time-dependent magnetohvdrodynamic (MHD) equations that describe 


( 1)187 
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atmospheric flows in three- dimensions for a spherical coordinate system can be written as follows: 


dp 1 d{r 2 pv r ) 1 d(pvesm0) 1 d(pv^) 
di r 2 dr r sin 0 80 r sin 0 8<j> 

d v r _ dvr_ _ vg dv T _ dv r _ 1 , d{pRT ) , dBe _ 1 dB r , 

dt dr r 80 t sin 0 d<f> p dr 9 dr r dO 

1 dB r _ dB+ v+ 2 + _ B^ 2 + B» 2 _ GM 

* r sin 0 d<$> dr + r pr r 2 

dv$ dv$ vq dv$ V+ dv$ 1 . 8{pRT) - , dBe l dB r 

dt Vr dr r 80 r sin0 d<i> p rd0 dr r 80 

1 8Be 8B+ . v r v* - v$ 2 cot 0 B r B$ - B$ 2 cot 0 
sin* 8<j> “ Td0 r pr 


dvj, _ dv+ vtidv^ dv+ B r dB+ Be dB* 1 

dt Vr dr r 80 r sin 0 d<f> + p dr r 80 rp sin 0 


dBr 

dt 

dB 0 

dt 


(tP + + B„^) + ^-(Br + Be cot 9) - — (t>, + v 9 cot 9) 

dtp dtp o<p pr r 

~ VaB ^ ~ T :( V + B r ~ v rB<t>)] 

r sm 0 80 dtp 

1 


7*T9lk {v9B * ~ v * Bs) ~ \b [r{VrB * - vaBr)] 


it m \b [r{v * Br ■ VtB * )] ~ l^e {v * B * ~ V * B * ] 

dT . 1 3(r 2 v r ) 1 d(?»smfl) 1 dv* 

dt ~ ^ r 2 dr + rsintf d9 rsin0 d<j> 

8T ve dT v+ dT 

Vr dr r 80 r sin 0 dtp 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


where the dependent variables are the density p, temperature 7\ velocity (u r , v$, v^), and magnetic field 
(B r , B<p). The independent variables are the radius r, the meridional angle 0 and azimuthal angle <j> as 
well at time ”t”. The constants are the polytropic index 7, solar total mass M and gravitational constant 
G. In addition, the standard equation of state (p = pRT) was used to obtain the above set of governing 
equations. 

The region, within which we will present the numerical solution to the above set of governing equations, is 
shown in Figure 1. This region is bounded by the solar surface and 15 solar radii ( R , ) in radial distance, 
by the equator and the pole in meridional distance (^-coordinate), and by azimuthal extent (^-coordinate) 
of 45°. 

Method of Solution 

The equations are solved numerically using a modified FICE (Full-Implicit-Continuous-Eulerian) scheme 
which is based on the original FICE scheme developed by Hu and Wu (1984); and Wu and Wang (1987). 
The grid spacings used are Sri — R,( 1 + M) 1 " 1 , and 60 — Stp = 4.5°. It should be noted that the radial 
spacing is not uniform and is chosen so as to; (1) assure the initial state as being in isothermal and 
hydrostatic equilibrium (Wang et al. 1982); (2) initialization of the computation procedure; and (3) to 
ensure numerical accuracy. The time step can be arbitrarily chosen because of the flexibility of the FICE 
scheme. 
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Fig. 1. A schematic description of the portion of a three- dimensional configuration in which the solution 
is calculated. Note that the computational domain extends from the pole to the equator within a 
45 ° extent of heliolongitude. In the present paper, symmetry is assumed below the solar equatorial 
plane. 


In itial State and Boundary Conditions 

In order to seek a solution of this problem, we need to specify the initial conditions These initial conditions 
include the magnetic field configuration, velocity field and corresponding thermodynamic properties of the 

plasma. 

The boundary conditions are rather complicated, hence a detailed account of the derivation of the boundary 
conditions will be presented later (Wang and Wu, 1990). We shall only briefly describe these boundary 
conditions here. There are a total of six sides in which the boundary conditions need to be specified, they 

q*)' r _ eight compatibility conditions are obtained from the set of governing equations (Wu and 
Wang, 1987); 

(2) r = 15fl,, non-reflecting boundary conditions are used (Hu and Wu, 1984); 

( 3 ) 0 = 0 (pole) and 9 = 90° (equator), symmetric conditions are chosen because of the chosen tiel 
configuration; 

( 4 ) <f> = 0 and d> = 45 °, the boundary conditions are obtained by extrapolation techniques. 

III. NUMERICAL RESULTS 

In order to carry out this simulation, we first introduced an initial state at isothermal and hydrostatic 
equilibrium with 7 = 1.67 together with a potential field in one case and, in a separate case, a linear torce- 
free magnetic field topology. These two separate cases were introduced into the set of governing e 9 uat >° n ® 
in order to ensure that the isothermal and hydrostatic equilibrium does exist. We then in ro 
steadv-state. Parker-type, velocity field. The numerical solution of this mathematical system led to a 
magnetohydrodynamic equilibrium state via the relaxation technique. This MHD eq num sta 
then taken as the simulated undisturbed coronal (i.e., quiet corona) with an outflowing solar win aroun 
multiple helmet magnetic topologies. 

The initial plasma and fields (magnetic and velocity) parameters incorporated in this simulation are the 
following representative conditions of a non-rotating sun with an initial plasma £>(= 16irn 0 jc„i», ) g 
unity, at r = R,, 9 = 90° and d> = 22.5 . 
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• Isothermal and hydrostatic equilibrium atmosphere. 

r„ = iO 6 K 

po = Po ex P> 1 7C ~ 1 


where pjj is the density (the value of 1.67 x 10 16 gm cm' 3 is used in this study and <?o is the {gravity on 
the solar surface. 

• Magnetic field configuration 

(i) A hexapole potential field (Jackson, 1962); and, in a separate calculation, 

(ii) A hexapole linear force-free field (Nakagawa et ad 1978) 

• Velocity Field 

v r (l, 0 , 6) — 15 km s _1 , 
u r (15, 0 , <i > ) = 200 km s” 1 , 
v$(ri 0, <p) = u*(r, 0, <j>) = 0. 

Figure 2 shows the simulated morphology of the quiet corona which consists of a three-dimensional repre- 
sentation of the brightness ( integrated density along the path of the line-of-sight), steady state solar wind 
velocity vectors and magnetic field for two cases: (a) initially potential field topology; and (b) initially linear 
force-free field topology, respectively. It is easy to recognize that the shape of the quiet corona depends on 
the initial magnetic field topology. The bright corona is related to the closed magnetic field configuration, 
and the dark region corresponds to the open field configuration which corresponds to the out- flowing solar 
wind from the coronal hole. Also it shows that the solar wind velocity is almost radial. 

In order to examine the physical structure of the quiet, steady-state, corona, we plot the radial distribution 
of the density and temperature at the pole and equator for the initially potential and linear force- free 
magnetic field topologies, respectively, as shown in Figure 3. The radial distribution of the three velocity 
components (i.e., t’n v^,) at the pole and equator is shown in Figure 4 for both types of magnetic field 
topology. Finally, we plot the radial distribution of Alfven and sonic speed at the pole and equator in 
Figure 5. 

Comparison of each of these parameters demonstrates the well-known inference and important fact that 
the magnetic field is the dominant factor that determines both the morphology and physical structure of 
the corona. The spatial diversity of these important, fundamental steady-state parameters is obvious. 

For the completeness of this presentation, we shall show some results for a disturbed corona in Figure 6. 
This numerical result is obtained by introducing a pressure pulse ( p/p 0 — 10) distributed over three grid 
points centered at 9 = 35°, 0 = 22.5° and r = R 9 for the case of the initially linear force-free magnetic field 
topology of the quiet corona as shown in Figure 2b. In Figure 6, at t = 600 s, we show simulated brightness 
(i.e. line-of-sight integrated density enhancement ), disturbed magnetic field and solar wind velocity vectors 
in the <p = 22° plane. According to the results shown, we may interpret that the brightness was caused by 
the flow interaction with the magnetic field. This density enhancement consists of both the mass carried 
by plasma flow motion and local wave compression. 

IV. CONCLUDING REMARKS 

In this study, we have presented a newly-developed, three-dimensional, time-dependent magnetohydrody- 
namic model for the study of corona structures in both quiet and disturbed states. This model extends 
from the solar surface to 15 R 9 and, thereby, includes the region of outflowing solar wind from the subsonic, 
sub-Alfvenic to super-sonic and super- ALfvenic regions. Therefore, we assert that it is, indeed, a model 
which could be used to study coronal/interplanetary coupling problems. 

In these preliminary results, we clearly recognize that the magnetic field topology and strength controls 
both the structures and physical parameters’ morphology of the corona. Also, this model has the capability 
to convert the fundamental physical parameters (i.e. p, T, t') to observables such as brightness (see Fig. 2) 
and doppler shifts (not shown). Therefore, we may claim that this model has the potential whereby it could 
be used as a diagnostic tool that can be applied to the interpretation and guidance of the observations. For 
example, we may use the physical properties obtained from this model to compute line profiles. As a final 
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Fig. 3. Radial distribution of the density (^/p Q )and temperature (T/T 0 ) at the pole and the equator, 
respectively for; (a) initially potential field configuration and (b) initially linear force- free field 
configuration with p a and T 0 given in page 6 of the text. 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

"s 

(a) 


Fig. 4. Radial distribution of the three velocity components ( v r , v$, v^) at the pole and equator for; (a) 
initially potential field configuration and (b) initially linear force-free field configuration. 
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Fig. 5. Alfvenic speed and sonic speed as function of radial distance at pole and equator for (a) initially 
potential field configuration and (b) intially linear force- free field configuration. 


remark, we recognize that the development of this model is far from complete. The improvements can be 
tackled in two major catagories as follows: 

0 Mathematical Improvement 

We should establish the accuracy of the numerical results. In order to achieve this purpose, we 
should conduct a grid size test for this model. 

• Physical Improvement 

Presently, the model includes dissipative mechanisms that were not invoked for the present demon- 
stration of its three-dimensional, temporal capability. Namely, the present model results are based 
on ’’ideal” MHD theory. We realize that dissipative MHD is important to many solar physics 
problems in which finite electrical conductivity, thermal conductivity, radiation and turbulence 
are undoubtedly present. We plan to incorporate these effects in our model via a conservative and 
rational step-by-step approach. However, the current ideal MHD model, because of its inherent 
and natural three-dimensional resemblance to the real world, is essential for the construction of 
solutions which resemble observed realistic topologies. We have obtained in the present demon- 
stration. for example, induced meridional and aximuthal flows which existing two-dimensional 
models cannot provide. 
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Abstract. Model calculations are presented for the rising motion of the top section oi a prominence loop, 
which is represented by a straight flux rope immersed in a coronal medium permeated with a bipolar 
magnetic held. Initially the prominence is at rest, in equilibrium with the surrounding coronal medium. When 
the maeneuc monopoies that account for the source current for the bipolar held strengthen, the upward 
hydromaenenc buoyancy force overcomes the downward gravitational force so that the prominence is 
initiated into rising motion. The illustrative examples show that prominences can move away from the solar 
surface by the action of the hvdromagnetic buoyancy force, which is preponderant with the diamagnetic 
force due to the current earned by the prominence interacting with the coronal magnetic field produced by 
the photospheric currents, if the changes in the photosphertc magnetic field are sufficiently large. 


1. Introduction 

A new dynamical model of prominence loops was recently constructed on the basis ol 
the theory of hvdromagnetic buoyancy force for fiux ropes (Yeh. 1989). A prominence 
loop immersed in the solar atmosphere is regarded as an extraneous body in the sense 
that it is magnetically separated from us surrounding medium. Thus us magnetic held, 
mass density, temperature, and motion arc quite different from those or the surrounding 
medium. The important feature is the polarization current induced on the periphery oi 
the prominence that makes the ambient magnetic field tangential. The exertion ot the 
ambient hvdromagnetic pressure gives rise to the hvdromagnetic buoyancy torce. Us 
predominant constituent is the diamagnetic force which amounts to the torce exerted 
on the currents in the prominence by the external currents that sustain the coronal 
magnetic field. For a prominence to be in stationary equilibrium with us surrounding 
medium, the hvdromaeneuc buoyancy force counterbalances the gravitational force 
exerted by the massive Sun. When the coronal magnetic field evolves, the changed 
diamagnetic force no longer matches the gravitational force. Once the torces become 
unbalanced, the prominence is initiated into motion. The evolving motion may be either 
upward or downward, depending on whether the hvdromagnetic buoyancy torce is 
greater or less than the gravitational force. That the evolving motion ot prominence 
filaments is driven by the evolution of the global magnetic field has been interred from 
observations (Kahier et aL 1988). 


Solar Physics 132: 335-351. 1991. 

D 1991 Kluwer Academic Publishers. Printed in Belgium. 
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la this paper we apply the dynamical theory to study the motion ot the top section 
of a prominence loop arched above the solar surface. The calculation presented is 
mathematically one-dimensional in space (viz., the heliocentric distance) although it 
involves two-dimensional geometry and the magnetic field is three-dimensional (see 
Figure 1). The governing equations for the dynamical evolution are MHD equauons or 
motion supplemented with equations of mass conservation, flux conservation, and 
energy conservation. Since we are mainly interested in the dynamics ot prominences, 
energetics is dealt with only to the extent necessary to provide a closed system ot 
equations for the dynamics. Accordingly, in our present calculations we neglect all 
entropy-generating processes that are pertinent to the thermodynamics ot prominences. 



Fig. I. Magnetic configuration resulting from the interaction between a couple ot magnetic monopoles on 
the solar surface and a large current earned by the prominence, with a polarization current induced on the 

interface. 


Several illustrative examples of dynamical evolution of a prominence loop are shown. 
First, we construct an equilibrium configuration for a prominence loop immersed in a 
coronal medium tnat has a bipolar magnetic field permeated into a magnetohydrostatic 
uimospnere. For a prominence loop to be stationary, it must have neither translational 
motion as a whole nor expansional motion relative to its axis. The former requires that 
the downward pull of the gravitational force exerted by the Sun is balanced by the. 
upward lift of the hydromagnetic buoyancy force exerted by the surrounding coronal 
medium. The latter requires that the outward push of the force by the gas pressure and 
the azimuthal current of the prominence is balanced by the inward pinch ot the force 
by the hydromagnetic pressure of the ambient medium and the axial current ot the 
prominence. Next, we calculate the motion of the prominence loop when the equilibrium 
is lost because of a temporal change of the bipolar field. The change can be caused by 
strengthening and/or displacement of the magnetic monopoles tor the bipolar field. 
These examples demonstrate that the prominence can move away trom the solar surtace 
when the change is sufficiently large. 

The model calculations illustrate the mechanism involving hydromagnetic buoyancy 
force that is likely to be important in the eruption of prominences. Such calculations 
in conjunction with analytical study also serve to narrow down the range ot the 
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parameters as an aid to MHD numerical simulations of the eruptive motion of promi- 
nences. Very often the difficulties with numerical simulations lie in the large number and 
extensive range of the pertinent parameters that characterize the phenomenon under 
study (Wu ei a/., 1990). 


2. Assumptions 

The geometry of the prominence loop may be described by its axis and its cross-section. 
We assume that the varying cross-secuon is well accounted for by a circular cross- 
section whose radius changes in time. In this treatment of the top section of a promi- 
nence loop, a prominence is represented by a flux rope with a straight axis, whose 
heliocentric distance may change. The corona is represented by a magnetized medium 
that has a transverse magnetic field, perpendicular to the axis of the prominence, which 
is bipolar and a longitudinal magnetic field, parallel to the axis, which vanes with the 
heliocentric distance only. The current that produces the bipolar field is below' the solar 
surface ; it is to be accounted for by a couple of magnetic monopoles on the photosphere. 
These monopoles are chosen to be line monopoles to make the problem two- 
dimensional. The current that produces the longitudinal magnetic field is in the corona. 
The coronal current is in magnetohydrostatic equilibrium with a stratified gas pressure 
of the coronal gas which is acted upon by solar gravity. 

The prominence loop cames helical field lines. The helical magnetic field in the 
straight prominence is represented by 

= 1 , 5 0 ^1 - — + ^<p 2^oQ ^ 

in cylindrical coordinates tr. q , <p), with the azimuthal angle cp measured trom the radial 
line pointing downward iu being the magnetic permeability in mks units). The axial 
component decreases from the axial value B n at the axis q - 0 to zero at the boundary 
q = Q. The azimuthal component increases from zero at the axis to the boundary value 
B b = kf-duQ at the boundary. The total axial flux is X ¥ E = i^Q 2 B {) and the total 
azimuthal flux is \Q 2 pJ n per unit axial length. This helical magnetic field is produced 
by the current density 


— l r /jj + 1 fpp. 


B» q!Q 

Q (I -q 2 iQ 2 ) 1 *' 


( 2 ) 


which has an axial component that is uniform and an azimuthal component that 
increases from zero at the axis to infinity at the boundary'. The total axial current is 
l E = nQ 2 J {) and the total azimuthal current is /i " 1 B 0 per unit axial length. The Lorentz- 
force density 


J E x - 1^ 



(3) 
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acting at various mass elements ol the prominence is in the radial direction, perpendicu- 
lar to the axis of the prominence loop. It increases from zero at the axis to 
pL-'B*iQ - ipJoQ at the boundary, in proportion to the radial distance. The axiai 
current produces a pinching force toward the axis whereas the azimuthal current 
produces an anti-pinching force away from the axis. 

The immersion of the prominence loop in the coronal medium causes a polarization 
current that keeps the internal field lines of the prominence separated from the external 
field lines of the corona. The induced current, which is concentrated in a thin peripheral 
layer bv virtue of the high electrical conductivity ot the solar plasma, produces a 
magnetic field that makes the ambient field tangential by cancelling the radial com- 
ponent of the coronal field on the interface and essentially doubling the azimuthal 
component there (Yeh. 1983). With the coronal mass density p.^, the coronal gas 
pressure p c , and the coronal magnetic field B , r s 1 . B . r . ~ B r x pre-existing at the sue 
of the prominence, the amoient magnetic field on the outer surtace oi the current layer 
is 


= l : B K „ mQ r 1*2(1*- B r .J, 


2kQ 


(4) 


The boundary magnetic field on the inner surtace is 


B B (<t>) = l 


^ E 

P 2kQ 


(5) 


The polarization current, given by i/, = l ? p“ l (B, - B 5 ) per unit circumterentiai length, 
shields off the coronal field from permeating into the prominence. Across the massless 
layer of peripheral current the sum of gas pressure and magnetic pressure is invariant. 
The ambient gas pressure 


pA<P) = Pr a 01 

on the outer surface is essentially equal to the pre-existing coronal gas pressure at the 
periphery' since the gas pressure in the exterior region is hardly perturbed by the intrusion 
of the prominence. The boundary gas pressure 


Pfl(0) “ P r. >q * fj 


X B\ 


5 CX q • Q COS" O 

+ <,-a cos * 

kQ 


(7) 


on the inner surface has a circumferential inhomogencity which is spatially transiormed 
from that of the ambient hydromagnetic pressure. The gas pressure inside the promi- 
nence is well represented by 

<t>) ~[Po - IP c + !«' ' B \ ).,-o] ~ -^ij + P rM ' + 


+ <P) + 4>) + 
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B, 


si _ q - c? 


COS <£ . 


U vanes from the axial value p (J at the axis to the boundary value p B ■ The gradient ot 
this gas pressure yields the iorce density 


-7 p E = -7 ( p x - ^-' B ^ + - \ r 


I E B ic. -L r - rn ^ 

* e 2 


2 Pn - <-P* + U_ ' B 2 r ) r-™ ( 9 ^ 

Q <2 ' 

ignoring lnsignincant terms. The term proportional to q;Q represents a radial torce 
density that resuits from the difference between the internal gas pressure and the 
circumferential average of the external hydromagneuc pressure. The otner terms, result- 
ine from the circumferential inhomogeneitv of the ambient hydromagnetic pressure, 
represent the spatial spreading of the hydromagnetic buoyancy force. The latter 
amounts to \ r p., GM 0 ir* + I r x B , r(2 2 + ' B\_ x ,R e by virtue of the magneto- 

hydrostatic state of the coronal medium (see Equations (24) and (25) for the definitions 
of R. and D- The coetficient f has the value of 2 in the above elucidation. 

The prominence moves with the velocity 


u f = U / 0 



( 10 ) 


which consists of a translational velocity common to ail mass elements ot the promi- 
nence loop and an expansional velocity proportional to the distance irom the axis. The 
translational motion is driven by the part of the force density that is unirorm and the 
expansional motion is driven by the part of the force density that is in various radial 
directions. The former part includes the gravitational torce exerted by the Sun and the 
hydromaenetic buoyancy force exerted by the surrounding medium. The latter part 
includes the Lorcntz force that results from the interaction among the internal currents 
inside the prominence and the gradient torce that results trom the pressure difference 
between the internal and external gases. 

The dynamical evolution of the prominence depends on its inertia. Vv e assume that 
the mass density is uniform over the cross section, ignoring the higher-order effect oi 
the spatial variation of the mass distribution. The value ot mass density p E changes in 
time. 


3. Governing Equations 

A prominence which is located initially equidistant Irom the two magnetic line 
monopoles will remain so when its heliocentric distance changes temporally (Figure 1). 
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The prominence loop is characterized by eight parameters: r 0 , Q , u 0 , V, p E , B n ,J n , ana 
P o • 

The charactenzing parameters evolve in accordance with the differential equations 


d 

— r 0 = m 0 • 
dr 




(ID 

i e =,: 

dr 




(12) 

d 

Pe 7" Ur > ~ Pi 
dr 

GM~ GM~ 

-J- , \ 4- 


R e 

(13) 

f: -» P r. a 

„ d v-'zEi 

*n"Br, 2p r - .a 

i 1 B 2 

- iuJ 2 Q 2 

1 14) 

dr 

Q 

0 



supplemented by the temporal invariances of total mass, axial magnetic riux. uzimutnai 
magnetic riux. and total thermal energy: 

nQ 2 p E = \(tr. V5) 

\«Q 2 B„ = 'fV , (16) 

iflV 0 -T/ e , (17) 

4 rc 


,)1 = 


(18) 


For a prominence to be initially tn stationary equilibrium with the surrounding coronal 
medium the requisite current density is 


J - B ‘ L ( | 9 ) 

7 Rr 

in terms of the mass density (or the requisite p E in terms of the current density) and 
other quantities. The requisite magnetic held, in either direction, is 

S 0 = - B\ - \prJlQr - 2pp n y 2 (20) 

in terms of the gas pressure (or the requisite p () in terms of the magnetic held) and other 
quantities. The first constraint makes the upward hvdromagnetic buoyancy force exactly 
balance the downward gravitational force. The second constraint makes the outward 
forces due to the gas pressure and the azimuthal current of the prominence exactly 
balance the inward forces due to the hvdromagnetic pressure of the ambient medium 
and the axial current of the prominence. These two constraints are depicted in Figures 
2 and 3, which show the required values of \pQJ 0 and 5 0 for various values oi r 0 > Q , 
P E , Po* P-n* TV, B X and B c m the neighborhood of the equilibrium values used 
in the examples (see Section o). 
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Fie. 2. Reauired values or' \uQJ„ and fl,, for an equilibrium prominence with various values ot r,.. ^ r . 

^>r p n . 


4. Coronal Medium 

We choose the line monopoles for the bipolar magnetic held to have strengths ot - *¥ sf 
and ancular separation of 2 U u . The two monopoles produce the bipolar magnetic held 

B = 1 B (21) 

n s. _ 

:n the midplane between them. Its direction, from the positive monopole to the negative 
monopole, is perpendicuiar to the midplane. Its magnitude is 

'j\r ^ 2 . sm —1 . (22) 

7T r 2 - 2rR^ cos t) u -r 

at a heliocentric distance of r. There the associated magnetic pressure has the gradient 






(23) 
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Fig. 3. Required vniucs ot UiQJ n ana 8 M for an eaumbnum prominence in a coronal medium uun various 

values oi ^ , . T , . 3 , . ur # , 


in the vertical direction, with the ratiius-of-eurvaturc 

^ r ~ -riR ^ cos -i- /? ~ p4) 

r - cos 

for the circular held line. (By virtue of the current-freeness of the bipolar held, the 
gradient force of its magnetic pressure is exactly opposite to its tensile force density.) 
This magnetic pressure gradient is enhanced by a factor 


u - £v'<n,)[(l - Q~'q\,r 4 *(? 2 'dJ 

by the polarization current 

i i -■ B zdMQl t 1 i1lLS 21 0 Z — i 

[Q 2 - 2Qq M Q0S[<p - <d u ) T q\ ( ] [ Q 2 - 2Qq M cositi + <b u ) t q\ f ] 

( 26 ) 
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Here = (r. ~ R c cos 0 Kf t Ri)' 2 is the distance from the prominence to either 
monopole and <p Kf = arc sin (R c sin 6 w iq Af ) is the azimuthal angle for the monopole. 
The field strength B _ x increases with and becomes maximized when cos d Kf is 
equal to 2 rR 0 j{r 2 4- R ^). In terms of the field strength on the soiar surface midway 
between the two monopoles, the monopoie strength has the value 


1 - cos 6*, 

Vu-2 * . ~ 

sin 


(27) 


The monopole strength [ F )f and the separation angle 20 > f may undergo temporal 
changes. 

We choose the longitudinal magnetic field 


B =1 : B 


(23) 


produced by the coronal current to be horizontal. The associated current density is 
I M |i~ 1 d B r dr. It provides a magnetic force in the force balance 


d , , , GMr 

" — iPr. - IP B-, ) « p, —r 

dr r- 


(29) 


between the gravitational force and the gradient of hydromagnctic pressure. In addition 
to the equation of force balance, two more constraints are needed in order to determine 
the vertical variation of the coronal mass density, gas pressure, and longitudinal 
magnetic field. We shall assume that the gas pressure varies in proportion to the mass 
density and the magnetic pressure vanes in proportion to the gas pressure, viz.. 


Pr. 


P T r ~ R- 

Pa ’ 

P r ' r « R - 


(30) 


P, -r-K 

These assumptions ensure that the pressure and mass density decrease with the helio- 
centric distance. They allow- us to calculate the mass density by numerical integration 
of the differential equation 


£ 1 GM~ d r 

1 r » Rr.( P r. • r » Rr* KT 7 ^ 


(32) 


trom the solar suriace. The ratio p f . ( divided by the gas constant of the solar plasma 
is the constant temperature T , of the coronal medium. 


5. Conditions for Upward and Outward Accelerations 

In order for the translational motion to have an upward acceleration away from the Sun. 
the hvdromagneiic buoyancy force must overcome the gravitational force. The former 
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will exceed the latter if the bipolar held is sufficiently large so that 


s. ri > I gjfs. 9 ' r - . (.33) 

i + li+ r5 7 0 

On the other hand, in order for the expansional motion to have an outward acceleration 
away from the axis, theoutward force must overcome the inward force. The former will 
exceed the latter if the surrounding medium has a hvdromagnetic pressure sufficiently 
small so that 


Px + " l5 >. 'I < Pn + zP ' B 2 - iP-JlQ 2 


(34) 


Upon the use of the equilibrium values ai / = 0 and the conservation invariants, the 
condition tor upward acceleration can be written as 



1 - 

l “ Px: -u 


(35) 


if we ignore the higher-order part of the diamagnetic force associated with the pre- 
existing gradient ot the coronal magnetic pressure. The condition for outward accelera- 
tion can be written as 


Px + zP~ 






(36) 


It follows trom the inequality ( 35) that in the reeion where p r is small, the translational 
motion will have an upward acceleration when the bipolar magnetic held encountered. 
B c x (r). is not less tnan its initial value by a factor of fr 0 t mih r 0 ) 2 . On the other hana. 
it follows from the inequality (36) that the expansional motion will have an outward 
acceleration in the region where p ^ ^ " 1 B l r is less than its initial value when 0{i ) 

is less than Q\ tm „ and in the region where p r r ] B 2 r . „ is sufficiently less than us 
initial value when 0(0 is greater than O 


6. Examples 

In mks units, the magnetic permeability has the value p = 4 k x 10 ” 6 7 T 2 m J J ~ *. the 
gravitational constant times solar mass has the value 

GM 0 = (6.67 x 10- M N m 2 kg- 2 ) x (1.99 x iO 30 kg), 

and the gas constant tor the proton-electron plasma has the value 

K = (1.38 x 10~ 23 J deg - l ):l(l-67 x 10 " 27 kg + 9.11 x 10' 31 kg). 
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To facilitate the numerics, we choose to measure time, length, and magnetic field in the 
units of 1 hour, one hehoradius, and 1 G. respectively, viz., f rcf = 3.6 x 10 3 s, 
• = 6.96 x 10 8 m, and B rc( = 10" 4 T. We further choose to measure speed, mass 


’ ref 


density, 
Cii 


pressure, current 


&rc f) ( r rcf/^rcf) 




in the units 


of 


(r re( lt Tef ) !K- Nameiv. 


z/ rcf = 193.3 km s 
irons cm" 3 ), 
T t ,r = 2.263 


- I 


densit> . 

£i ' 

2.129 x 10' 

3 ), /? rcf = 7.958 x 10' 3 J m -3 , J rcf = 1.143 x 10” 7 A m 2 , and 
x 10 b deg. in these normalized units, henceforth, both the magnetic perme- 
ability. aiven bv 


and lemperature 
‘ B r Jr rel . and 
13 kg m " 3 (corresponding to 1.274 x 10 3 elec- 


(4 it x 10" 7 T 2 m 3 J 


^ref/^ref^rer') 


and the gas constant, given by 

(1.65 x 10 4 J kg “ 1 deg ' 1 ) (u 2 cf , T rct ), 

have the value of unity whereas the gravitational constant times soiar mass has the value 
of 5.102 helioradius 3 h " given by (1.333 x 10 2o m 3 s “ 2 )/(r ref w 2 ef ). 

For the magnetohvdrostatic coronal atmosphere, we choose a mass density of 
3 x 10 9 electrons cm ~ \ a temperature of 2 x 10° K (hence, 7~ x = 0.8838) and a longi- 
tudinal magnetic field of 2 G at the solar surface. For the magnetic monopoies. to have 
a transverse magnetic field of 10 G at the soiar surface midway between the two 


gams 



Fig. 4. Profile of a stratified magnetohvdrostauc corona at 2 x 10° K. with a mass density of 
3 x 10 V electrons cm \ a longitudinal magnetic field of 2 G, and a transverse magnetic field of 10 G at the 

solar surface. 
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monopoles, we choose 

= 5.497 G helioradius. 0 Kf = 10 3 . 

The calculated prohies are shown in Figure 4. 

For a stationary prominence loop, we choose a height of 5 x L0 4 km, a radius of 
2xl0 4 km, a mass density of 5 x 10 11 electrons cm “ 3 , and a temperature of 
5 x 10 4 K so that: 

r 0 - 1.0718. Q = 0.02874. = 3924.2, p 0 = 86.704. 

At the site where the prominence resides we have 

= 16.535. p r = 14.614. 3 r , = 1.676. = 8.054. 

The conditions of force balance require 
J t) = 2080.9, 3 0 = 40.574. 

In other words, for the prominence loop to be in stationary equilibrium with the 
surrounding medium, it must carry a total axial current / E of 3.0 x 10 M A and carry an 
azimuthal current that sustains a total axial magnetic dux K ¥ E of 3.4 x IQ 12 webers. 
These values are within the range of typical values for quiescent prominences 

ip 


A/m 



Fig. 5 . Peripheral distribution 01' the polarization current in the equilibrium contiguration. 
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(Tandberg-Hanssen. 1974). It is seen from 


Pe 


GA/ 2 
>*o 


17428.7. 


GM, 


73.61, 


1eB ^ ± = 16759.7, 

*Q 2 


R c 

that the gravitational force is largely balanced by the diamagnetic force (the part due 
to the inhomogeneuv of the coronal magnetic field is only 3.55° 0 of the part due to the 
prominence current ). The hydrostatic buoyancy force is very small, accounting ior only 
0.42° n . On the other hand, it is seen from 

P{) = 86.704. iu-'S-l = 823.14. \pJlQ 2 = S93.82. 

p, = 14.615. = 1.405 



helioradu /hours he hondu/ hours 




Fig. 6. Evolution ot the prominence as the monopole strength increases, 
with d¥\,;di - 20 G helioradius h" 
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that the pinching force of the axiai current is largely balanced by the anu-pinchmg force 
of the azimuthal current and to a lesser extent by the internal gas pressure. The ambient 
hydromagneuc pressure provides only a very small pinching. To facilitate comparison, 
these values may be translated to (2/ip x / 2 = 5,4063 G, ~ -9.903 G. and 

(2#ip 0 ) l;2 = 13.168 Q. The plasma beta at the axis is 0.1053. The transverse projection 
of the field lines in the equilibrium configuration is as shown in Figure 1 . With the ratio 
M/j/lPi, = 0.9684, the bipolar field has two neutral points located at q = 0.187. 
<p = ± 1 1 8.7 c . Figure 5 shows the peripheral distribution of the polarization current with 
q M = 0.194 and <p s1 = 63.4°. The polarization current flows in the direction of the 
prominence current in the lower periphery ' <p\ < 82.6~ and flows in the opposite 
direction in the upper periphery. It is zero at the two points where the two neutral points 
would be located in the case / F happens to be zero. Of course, the total polarization 
current sums up to zero. 


hebondu 



hcLiondu/houn 



Fig. 7. Evolution ol the prominence with various values ot df„/d l. 
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Now, we consider the motion of the prominence loop when the equilibrium is lost. 
The initialed motion is caused by a temporal change of the bipolar magnetic field due 
to the strengthening of the monopoies. Figure 6 shows the evolution with 

— *P Xf = 20.0 G helioradius h " ] for 0 < i < 10 h , 
dr 

revealed by the numerical solution of the dynamical equations. The initial increase in 
the bipolar field makes the hvdromagnetic buoyancy force exceed the gravitational force, 
so that the prominence rises from its equilibrium position. The prominence keeps 
moving upward, even'during the short while 0.6 < r < 3.5 when the prominence is 
decelerated because the hvdromagnetic buoyancy force is not large enough to exceed 
the gravitational force. Likewise, the radius of the prominence keeps increasing. Its rate 
of increase is small in this case because the encountered coronal hvdromagnetic pres- 
sure decreases very siowlv. To see the dependence on the speed of the evolution of the 
photosphenc magnetic held, we show in Figure 7 the evolution caused by smaller 
values of d i. It is seen that the prominence may move up and down in oscillational 
motion if the evolution is slow. With a sufficiently rapid evolution, the prominence will 
move away from the Sun. In passing, we mention that the prominence would move 
downward in response to negative values of 


7. Discussion 


The calculations show the importance of the hvdromagnetic buoyancy force, which 
includes the diamagnetic force, in the dynamics of prominence loops. In the illustrative 
examples, not only in equilibrium but also during the motion, the hvdromagnetic 
buoyancy force is preponderant with the diamagnetic force due to the current earned 
by the prominence interacting with the coronal magnetic field produced by the photo- 
spheric currents. The part of the diamagnetic force due to the inhomogeneity of the 
coronal magnetic field only amounts to a few percent and the hydrostatic buoyancy force 
is even much smaller. Their percentages diminish in hcliocentnc distance. Therefore, we 
conclude that without the action of the diamagnetic force due to the prominence current, 
prominences are not likely to move aw'av from the solar surface. 

In fact, when the prominence is close to the photosphenc surface, the diamagnetic 
force is enhanced by the mirror-current effect (Kuperus and Raadu. 1974). The polari- 
zation current induced on the photosphere will exert an additional upward force on the 
prominence current in the amount of fiI 2 £ R nr 0 lr^ - R%) (cf. Van Tend and Kuperus. 
1978). Inclusion of this force will modify Equation (13) to 


d 

% 7 - f 'o 
dl 


-Pt' 


GMr- 


GM, 


JnB^ ± - r 


M-'S- 


R, 


+ nJfi 


Q^Rl 

- r<Arf, - Rk) 


(37) 
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Fig. 8. Evolution of the prominence when the mirror-current effect is included in the hvdromaenetic 
diamagnetic force, with d ¥\ r /d/ ^ 20 G heiioradtus h 1 . Dashed lines indicate the corresponding evolution 

without the mirror currents. 


Accordingly, for the prominence in the illustrative examples, the requisite current density 
reduces to J 0 = 1426.5 and the requisite magnetic field reduces to B n = 26.44. With the 
monopoles strengthening at the rate of dV Uf f dt = 20. the prominence rises slightly 
slower. See Figures 8. This is due to a smaller prominence current. The reduced l E B x x 
is not sufficiently compensated by the added sr 0 (rg - R 2 0 ). This calculation 

serves as an a posteriori justification for the simplified model by Van Tend and Kuperus 
(1978). 

The dynamical model illustrated by the calculations presented in this paper applies 
to prominence eruptions driven by gradual evolution of the photospheric magnetic field 
in weak field regions. It also applies to filament eruptions in active regions, with other 
parameter values. 
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